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ABSTRACT 

A conception of virtual quantum information bit - virtual qubit - is introduced. It 
is shown by means of virtual qubit representation that four states of a single quantum 
particle is enough for implementation of full set of the gates, which is necessary for 
creation an arbitrary algorithm for a quantum computer. 

The physical nature and mutual disposition of four working states is of no 
significance, if there are suitable selection rules for the particle interaction with the 
external electromagnetic field pulses. 

1. INTRODUCTION 

The search of new physical systems, capable to be a medium for implementation 
of the quantum computation, is one of the direction which promises appreciable 
progress in quantum computer realization. For this aim it is necessary to find an 
adequate quantum system and assign qubits to a subset of its stationary states. As the 
next step one should find the physical means which cause the transitions between these 
states, realizing quantum gates. Besides that it is necessary to determine how to create 
initial state and how to readout final state. 

Spin 1/2 nuclear magnetic resonance is one of the examples of such physical 
systems 1 ' 2 . Two possible spin 1/2 stationary states represent naturally one information 
bit (quantum bit or qubit). Thus generally used quantum computer model has the 
following correspondence: one qubit - one two level particle. The possibility to use the 
arbitrary physical origin quantum object possessing four or more discrete energy levels 
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as a basis element for storing and processing information in a quantum computer is 
considered in the present paper. Main results have been published 3 . 

In standard quantum computer models the exchange interaction is used to realize 
two qubit gates. To implement the gates it is necessary to turn on the particle interaction 
for accurate time intervals and to turn it off when a gate should not operate. These time 
intervals are determined by interaction constant and may be longer than coherence time. 
In addition, to turn the interactions off the application of the complicated experimental 
schemes are necessary which can make the computation time very long. 

In proposed model the gates are implemented on the single quantum particle and 
it is not necessary to spend time for the distant particle interactions have worked out. 
The logic gates in the proposed model are realized using external field short pulses. The 
pulse duration are defined by the field amplitude and are under full experimental 
control. 

2. PHYSICAL SYSTEM 

To be concrete let us consider four level discrete spectrum of a nucleus with spin 
3/2 and quadruple moment, which is placed in a constant magnetic field and crystal 
electric field gradient. 

2.1. The main Hamiltonian 

The spin energy levels are defined by the Hamiltonian: 

Ho = H z + H Q , (1) 

where H z - Hamiltonian of interaction with external constant magnetic field , H Q - 

Hamiltonian of electric quadrupole interactions. In particular case, when constant 
magnetic field is parallel to one of the electric field gradient principal axes, these 
operators take the form 4 
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H z = - too / z , H Q = l*co Q [3/ 2 - /(/+!)+ ri (7 X 2 - / y 2 )] (2) 



where Cfl Q - quadrupole interaction constant, (Oq - Zeeman frequency, T] - the electric 

field gradient asymmetry parameter (lr|l<l). In the spin 3/2 case the Hamiltonian (1) ,(2) 
eigenvalues are equal to 



h 8+3/2 = h to Q (B T + C) ; ft e±i/2 = h <o Q (- B± ± C); (3) 



where 



B+ = V [(1±2C) 2 + r| 2 /3]; C = co Q /(0 Q ; 
The corresponding eigenfuctions are 

3/2 >= C0s(CX±) |% + 3/2 >+ 5m((X+) |%±l/2>; 

ra*± 1/2 >= co5(a±) |% ±1/2 > - sin(a±) |% T 3 /2>; (4) 
where l% m > is the operator I z eigenfunction corresponding to eigenvalue m and 
fgctt=(V3) [B± + (l±2C)]/ri. 

To simplify the notation the indexes -3/2, -1/2, +1/2, +3/2 will be replaced by 1, 2, 3, 4 
correspondingly. The above formulas (l)-(4) are valid at arbitrary ratio R = G\/(Qq- To be 
concrete the case when the magnetic resonance spectrum is split by quadrupole 
interaction (co >Cfl Q ), at the condition that G) Q is much grater than the spin level width is 
considered below. In this case the spin spectrum consists of the four well resolved 
resonance lines with energies £i>£2>£3>£4. 
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2.2. The resonance radio frequency field interaction Hamiltonian 

Neglecting the relaxation, the full nuclear spin Hamiltonian can be written as 



H = Ho + H rf , (5) 

where H rt is the operator of interaction with RF field. In the case when the RF field 
vector is parallel to Y axis this operator is 



H Y d (t) = 2fiy Hrf/ Y Cos(Qt +/), (6) 



where Hrf , / h Q. are the RF field amplitude, phase and frequency, y is the nucleus gyro 
magnetic ratio. 

It will be needed later the operator HJt) interaction representation form: 
H*Jt) = D- l (t-t ) H rf (t) D(t-to), 

D(t-to) = exp [ - /H (t - to)/* ], (7) 

The transformation operator D form corresponds to the case when the main Hamiltonian 
is time independent. 

The representation of projective operators, P mn , for spins will be used to simplify 
notations and calculations. Each of these operators P mn is 4x4 matrix with all zero 
elements except one p mn = 1. Here m,n = 1,2,3,4 and thus one has the set of 16 
operators. The spin components are expressed in the projective operator terms as 

<¥ m IJ„ I^Pmn 

The projective operators have very simple multiplication rules 

PklPmn = Sl m Pkn 5 Pmn = P nm 
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Pmn| l i / k> = 5 nk |^ m > (8) 

In the projective operator representation the main Hamiltonian has the form 

H„ = E m M m P mm , (9) 

and the transformation operator is 

D(t-to) = E m P mm exp [ - /em (t - 1 ) ] . (10) 
The RF interaction Hamiltonian can be represented as 

H* rf (t) = H d eff + Z kl G kl exp[ /(8 k -8i )t] (11) 

when the RF frequency fi coincides with one of the spin eigenfrequencies, for example 
with Q mn = le ra - 8 n l. All operators in the equation (11) are of the same order of value: 

|G k il ~ |H rfj e ff | ~ |H rf |, and both G k i and H rf> eff are time independent. The last terms in 
(11), having amplitudes |G k il, oscillate with eigenfrequencies Q k i ■£ Q. mn . Their influence 
on the spin evolution can be neglected in comparison with that of the constant term 
|Hrf ; eff |. This is due to the fact that in NMR experiments the RF pulse duration is usually 
much more than the cycle period. Thus the interaction (6) excites effectively resonance 
transitions between the energy levels with Q. = Q. mn only. That is why the interaction 
representation operator (6) is reduced over several cycles (Q k i t »1) to effective 
operator 

H Y rf eff = i tq H rf k¥ m l/ Y l¥ n >l(P mn j f - P nm Q if ). (12) 

It is necessary to point out that although in the basis l%> the operators / , / have 

matrix element selection rules m-n = ±1, in the basis \ x ¥> these operators have non zero 
matrix elements with 
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m-n = ±l,±2. (13) 

2.3. The spin system evolution under RF pulses 

The state vector evolution 

W(t)> = U(t,to)W(to)>, (14) 
of a physical system can be defined using unitary evolution operator 5 : 

U(t,t ) = D(t-to)V(t,t ), (15) 

V(t,t ) = T exp [ - (M) J H* rf (f) dt'], 
H* rf (t) = D" 1 (t-to) H rf (t) D(t-to), 

where T is chronological Dyson operator and V is the evolution operator in interaction 
representation. 

At the condition, when the fast oscillating terms in H*rf(t) can be neglected, the 
chronological exponent transforms to ordinary one and the operator V(t,t ) reduces to 

t 

V Y (t,t ) - exp[ -ilfi J H Y rf , eff dt'] = exp [((p/2)(P mn e lf - P nm e" lf )] = V Y (a mn , (p y ). 
where cp y = 2(t-t )YH rf |< K F n l/ y l*P m >| and e m > 8 n . 

Expanding the exponent and using the projective operator multiplication rules (8) one 
has 



V Y (a mn , q>y) = P kk + P u +(P nn + P rara )co5((p y /2)+ (P mn e if - P nm e" if ) sin(%/2). (16) 



Here the indexes k,l ^ m,n. The evolution operator expression Vx(^ mn , 9x) for the case 
when RF field is parallel to X axis can be obtained from Eq.(16) substituting the phase/ 
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by (f- nil). In the present paper all necessary quantum gates can be realized with phase 

f=o. 



To implement quantum gates two frequency excitation will be necessary too. 
Below it will be shown that excitation of two transitions which have no common energy 
levels are enough for this aim. In these cases the evolution operator is a product of two 
commuting evolution operators which belongs to different transitions. For two 
excitation frequencies - Q. n and Q 34 - the straightforward calculations give 

V Y (Q 12 , cp Y ; ^34, 9'y) = (Pn + P22) cos(<p Y /2) + (P 21 - P 12 ) sm(<p Y /2) + 

+(P 33 + P44) cos(<?' Y /2) + (P 43 - P 34 ) sin(<?' Y /2). (17) 

For two excitation frequencies - and Q. 24 - the evolution operator is 

Vy(£2i3, (p Y ; ^24, 9'y) = (P22 + P44) cos((p Y /2) + (P 42 - P 24 ) sin((p Y /2)+ 

+(P 33 + P11 ) cos(q>' Y /2) + (P 3 i - P13) sin{a?' Y l2). (18) 



3. THE VIRTUAL SPIN FORMALISM 



1 2 

In the adopted NMR computer standard model ' the basis for quantum gates 
implementation consists of two real spins R=l/2 and S=l/2 with exchange interaction. 
The states of such a system are written usually in the quantum theory language as the 
abstract four dimensional space which is the direct product TR®rs of real spin two 
dimensional eigenstate spaces T R and T s . To make clear the information aspect of 
suggested logic operations in our case it is useful to apply the inverse procedure: to 
represent real spin 3/2 four dimensional space Y\ as a direct product r R ®r s of two 
abstract two dimensional spaces of the virtual spin R and S states. Then any four 
dimensional operator P can be expressed as a linear combination of direct product R®S 
components. There is the following isomorphic correspondence between the space Ti 
basis I v I / m> and the direct product r R ®r s basis l^ m >®l^ n > 
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Wl> = l£l>®l(;l>slll>, 

I 1 I / 2>=I^1>®IC2> = I10>, 



l^ 3 > = l^2>®lCl> = I01>, 
l X P 4 > = I^2>®IC2> = I00>, 



(19) 



where virtual spin indexes -1/2 and +1/2 are replaced with 1 and 2, correspondingly. 
Here II 1> and so on are the notations, adopted in quantum information theory to 
represent the two qubit states. There are the following useful relations between the 
different spaces projective operators 



1r ® S k lS mn = 8l m l R ® S kn ; 1r ®S k l l^ ra >lCn> = 5l n l^ m >lCk> ', 

(Rki®S a b)( R mn ®S C d )= 8i m 8bc Rkn ® Sad ; 

Here the operators R mn and S mn are the Tr and Ts space projective operators whereas 1r 
and Is are the unit operators in the corresponding spaces. In particular the virtual spin 
components R x , R y , R z and S x , S y , S z are expressed as 

R x = (R 12 +R 21 )/2; R y = i(R 12 -R 21 )/2; R z =(R u -R 22 )/2; 
S x = (S 12 +S 21 )/2; S y = /(Si2-S 2 i)/2; S z =(Sn-S 22 )/2. 

It is necessary to point out the fact that in four dimensional space Ti two qubits can be 
defined in another way. One can consider that the upper two levels compose the first 
qubit whereas the bottom pair of levels composes the second qubit. Such two level 
system definition is commonly used in magnetic resonance under the name ^fictitious 
(or effective) spin 1/2 formalism*. In this case the space Ti is a direct sum Tr©!^ of 
fictitious spin spaces. Using such an approach one can easily implement single qubit 
rotations, but the two qubit controlled NOT» gate implementation will face with some 
problems. Besides that the introduced virtual spins correspond very well to quantum 
information language. 




(20) 
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4. INITIAL STATE PREPARATION 



The state I00> is the initial input state for present days quantum algorithms. To 
emulate an abstract quantum computer by means of NMR experiments it is necessary to 
take into account the following peculiarities of NMR quantum computing. In NMR 
information processing the macroscopic number of identical quantum processors - 
molecules are used in parallel and the output is a sum of these molecules' signals. The 
density matrix formalism is the adequate language for NMR experiments. The density 
matrix G imt = const P44 is the equivalent of the I00> state. 

The macroscopic sample of identical nuclei ensemble has energy levels with 
population distributed according to Boltzman law 

p eq = Z" 1 exp( - p H ), Z = Sp[exp( -$H)], (3 = I/kT, (21) 

In the usual NMR experiments lp\ff l«l (normally ~10" 5 - 10" 6 ). At such conditions the 
input state for NMR quantum computer is given by matrix 

Peq = Z [lj + 52 m A m Pmm] 

li = £m Pmm; A m = fe m /kT; (22) 

It is possible to prepare p eq in the form Gmn directly by cooling spin system till ultra low 
temperature. It will create not only tremendous technological difficulties but also will 
slow down the computation speed. Below it will be shown how in real NMR 
experiments to get the density matrix 

pinit = const [li + const P44] (23) 

as an input for quantum algorithms. In comparison with G init there is an additional term 
in pi n i t which is proportional to unity matrix li .But it does not contribute to observable 
signal being invariant under unitary transformations. Therefore any pulse sequences, 
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including the computational ones, do not affect the matrix li and that is why the initial 
matrixes p^t and Oinit give the same output result. 



The space 2 or temporal 6 averaging procedure can be used to transform the 
density matrix p eq to the desirable form p init . Let us to consider temporal averaging 
procedure for multiqubit spin, leaving other methods for future developments. 

The following procedure, which goes back to the paper 6 ideas, can be suggested. 
Let the desirable computation consists of the fulfilling the transformation V com p with the 
initial state Pi mt , whereas spin system has an equilibrium density matrix of the form (22). 
It can be shown that the average 

(l/3)[V com pP e q V comp V com pVi p e q V lV c omp V com pV2 Peq V 2 V C omp ] — 
=(l/3)V comp ( p e q + Vl peq V\ + V 2 peq V + 2 ) V + 

comp 

of three computations with the same propagator V comp , all of which start from three 
different input states p eq , Vip eq V + i and V2p eq V + 2 , is equivalent to transformation V com p 
of the state p init : 

VcompPinit V c omp 

(24) 

To prove this statement we choose the transformations Vi and V2 in the form of two 
successive single frequency pulse propagators 

Vi = v Y (a 12 ,7i)VY(a23,7i) = p 44 + P21 + Po + P32, 

v 2 = v Y (a 2 3,7i) v Y (a 12,71)= p 44 - P12 + P31 - P 23 (25) 

Using equations (16a), (16b) and multiplication rules (8) we get 

(1/3) [p eq + Vl Peq V + l + V 2 Peq V + 2 ] = Z [0C 1 + (3P 44 ] - Pinit , (26) 
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a = l+(l/3)[Xi +l 2 + hl (3 = X 4 - (1/3) [Xi +l 2 + hi 



This proves the possibility to start with p eq instead of pMt using suitable pulse 
sequences. 

5. ONE QUBIT ROTATION 

The operator Vy(^i2, cp Y ; ^34, <p'y), Eq. (17), at the condition cp Y = <p'y = cp is 
equal to 

V Y (Q 12 , (p; Q 3 4, 9) =(Pll+P22+P33+P44)C05((p/2)+(P 21 -P 12 +P 4 3-P34)«'«(Cp/2). (27) 

Expressing the operators P in terms of operators R and S, Eq. (20), it can be shown that 
the operator 

V Y (Q 12 , cp; Q 34 , cp) = (Rn + R2 2 )®[ (Sn + S 22 ) cos(<?/2) + (S21 - Si 2 )sm(cp/2)] = 

= exp{ - i (p l R ®S y } (28) 

is equivalent to operator (27). The last equality in the right hand side of (28) means that 
the transformation Vy(^i2, cp; ^34, cp) is a rotation through the angle cp around Y axis in 
the space T s keeping the space T R invariant. 

Using the same arguments it can be shown that the operator V Y (^i3, cp Y ; ^24, cp' Y ), Eq. 
(18), at cpY=cp'Y=cp is equal to operator 

V Y (Q 13 , cp; ^24, cp)= [(Rii+R22)co^((p/2)+(R 2 i-R 12 )5m((p/2)]®(Sn+S 22 )= 

= exp{ -i(pR y ®l s }. (29) 

The last equality in the right hand side of (29) means that the transformation V Y (Qi 3 , cp; 
Q.24, cp) is a rotation through the angle cp around Y axis in the space T R keeping the space 
T s invariant. 
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6. CONTROLLED NOT» GATE - CNOT 



The transformation V y (£2i 2 , cp y ), Eq. (16), which corresponds to energy levels h£\ 
and he 2 and is defined as 

V y (Q 12 , 71) = [P 33 + P 44 ] + [P21 - P12], (30) 

at cp y = 71, implements two qubits operation controlled NOT» - CNOT: NOT operation 
at spin S if R is in state l^i> and leaves S in its current state if R is in state l^ 2 >. 
Straightforward calculation using (30) gives 

V y (Q 12 , 7C) l^!> = ^ V y (Q 12 , 7C) l¥ 2 > = - FJ?i>, 

Y y (Qn, 71) l v I , 3> = l v i , 3>, V y (Q 12 , 71) W 4 > = l^4>, 

or in information notations 

V y (Q 12 , 71) II 1> = I10>, V y (^i 2 , 71) I10> = - II 1>, 

V y (Q 12 , 71) I01> = I01>, V y (n 12 , 71) I00> = I00>. 

It is the truth table for CNOT operation defined above neglecting phase factor exp(i%). 
The evolution operator V y (£2i 2 , n) in the basis T R ®rs projective operators representation 
can be written as 

V y (Q 12 , 7C) = R„ ® (S21 - S12) + R22 ®ls (31) 

Using the same arguments it can be shown that transition transformation V y (£2i 3 , cp y ) 
between h 81 and ^283 states at cp y = 7t 

V y (n 13 , 71) = [P 22 + P 44 ] + [P31 - P13], (32) 
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implements CNOT operation in which the spin S state controls the spin R state. 
Operator (32) gives the truth table 



V y (Q 13 , 71) l¥i> = l^3>, V y (Q 13 , 7C) l¥ 2 > = l^ 2 >, 

V y (Q 13 , 71) W 3 > = " l^l>, V y (Q 13 , 71) l^ 4 > = W*>. 

which in the information notations has the form 

V y (£2i 3 ,7i)lll> = I01>, Vy(fii 3 , 7t) I10> = I10>, 

V y (Hi 3 , ji) I01> = - II 1>, V y (Q 13 , 7t) I00> = I00>. 

Using the basis r R ®r s projective operators the evolution operator V y (£2i 3 , n) can be 
written as 

V y (Q 13 , 7t) = 1 R ® S 22 + ( R21-R12) ® Sn. (33) 

Thus the spin 1=3/2 pulse excitation corresponding to the evolution operators V y (ni2,7t) 
and V y (Qi 3 ,7C) implements the quantum CNOT gates of two kinds - when spin R 
controls the spin S dynamics and vice versa. 

7. READOUT 

To obtain results of computation it is necessary to have possibility to read out the 
density matrix final state p out . The NMR methods allows to get all matrix elements using 
so called «state tomography» 7 . This method is realized by means of complicated pulse 
sequences. For other physical systems the development of adequate means is necessary. 

8. CONCLUDING REMARKS 



13 



The presented formalism is applicable to quantum system of arbitrary physical 
nature. One can use the four appropriate energy levels of greater nuclear spins, electron 
paramagnetic resonance spectra with effective spin S* > 3/2, optical levels and so on. 
Only the resonance frequency values and the operators matrix elements will be changed 
in different physical realizations. For example, the nuclear spin 1=3/2 quadrupole 
resonance spectrum being split by interaction with constant magnetic field (co < G) Q ) can 

be used. This case differs from the one considered in the paper only in mutual energy 
levels arrangement whereas the equations (1)-(10) for operators and eigenfunctions 
remain the same. 

The virtual spins representation provides a new way to implement three qubit 
gates. In future such gates can be useful for development the compact algorithms or for 
other purposes. But its implementation using two level quantum systems requires three 
body interactions which are absent in nature. One can notice that spin 7/2 Hilbert space 
Ti can be considered as three virtual spin 1/2 Hilbert spaces direct product. Therefore the 
Hilbert space of one such particle can contain three qubits. Then, finding suitable pulse 
sequences that act onto eight dimensional spin 7/2 system, it is possible to implement 
the gates, for realization of which three individual spins 1/2 interaction is necessary. In 
another solution of this problem was found and it was shown that special two qubit gate 
set based on two body physical interactions is enough to realized any arbitrary quantum 
algorithm. 
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